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Abstract: 

The statistics of S-matrix fluctuations are numerically investigated on a model for 
irregular quantum scattering in which a classical chaotic diffusion takes place within 
the interaction region. Agreement with various random-matrix theoretic predictions is 
discussed in the various regimes ( ballistic, diffusive, localized ) 
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Erratic fluctuations of scattering cross sections in dependence of energy or other con- 
trol parameters are a well known phenomenon in several branches of quantum physics. 
Stochastic methods which model the Scattering Matrix by a random matrix chosen in 
an appropriate ensemble have proven very successful in analyzing the statistics of these 
fluctuations. These methods, originally started in Statistical Nuclear Physics, are now 
widespread in other fields of microphysics[l]; of particular importance among these, the 
analysis of conductance fluctuations of small metallic samples at low temperature. A pos- 
sible unifying approach to such disparate fluctuation phenomena has emerged from studies 
on the quantum dynamics of classically chaotic systems. Indeed, recent developments in 
that field indicate that cross section fluctuations distinguished by some universal features 
generically appear in quantum scattering processes when the corresponding classical dy- 
namics is chaotic inside the interaction region[2, 3]. In some model systems exhibiting 
classical chaotic scattering, the statistics of the rapid S-matrix fluctuations occurring on 
energy scales comparable to the mean level spacing have been found to closely agree with 
the predictions of Random Matrix Theory, and specifically with those of Dyson's Circular 
Orthogonal Ensemble (COE) in the case of time reversal invariant scattering processes. In 
other cases, more specific matrix ensembles have been constructed in order to account for 
the presence of fast direct processes besides the chaotic one, and again a good agreement 
has been found [4, 5]. These results suggest that classical chaos may eventually provide 
dynamical grounds for the use of Random Matrix theoretical methods. 

Analyzing conductance fluctuations from the standpoint of quantum chaotic scattering 
is a promising task. Some significant results in this direction have been obtained for the 
ballistic regime[6], but to the best of our knowledge the fluctuations of the S-matrix 
have not yet been investigated in situations exhibiting at once classical chaotic diffusion 
and quantum localization. This is the subject of the present Letter, where this issue 
is investigated on a recently introduced one dimensional model that can be analyzed 
both classically and quantum mechanically [7, 8]. The main qualitative features of this 
model are as follows. In the classical case it displays a chaotic diffusion leading to ohmic 
conductance. In the quantum case, the interplay of diffusion and localization determines 
a transition from a ballistic regime, where the sample length is comparable to the mean 
free path, to a localized regime, in which the sample length is larger than the localization 
length. At intermediate lengths between these two regimes, the quantum transmission 
coefficient exhibits an ohmic dependence on the sample length. This crossover region 
increases in size on moving towards the classical limit, and within it the magnitude of 
transmission fluctuations is scarcely dependent on the sample length. Thus the model 
exhibits some essential features of classical and quantum transport, and the analysis of its 
scattering fluctuations is relevant for the above sketched general themes. In this Letter 
we compare the numerically computed statistics of the S-matrix with the predictions of 
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random matrix ensembles. 

We shall now briefly describe the model, which was presented in detail in refs.[7, 8]. It 
describes the discrete time ( stroboscopic ) dynamics of a particle on a line. The scattering 
region ( in solid state language, the "sample" ) is an interval of length L. The one-step 
classical evolution is defined by: 

n = n-V'{9) 9 = 9 + rn (mod(2ir)). (1) 

where n defines the position of the particle on the line, V(9) = k cos 9 with k a fixed 
parameter, r is zero outside the scattering region defined by M : no < n < L + no an d 
has a constant nonzero value ( also denoted by r ) inside the region itself. The interacting 
dynamics is thus defined by the Standard Map; the free dynamics is just uniform motion 
at a constant speed k sin 9. Changes of the phase 9 occur only inside the sample and can 
be thought of as "collisions". The mean distance travelled between subsequent collisions 
is on the order of k which therefore defines the mean free path (mfp). In the deeply 
chaotic regime of the Standard Map (kr >> 1) chaotic transport takes place inside the 
sample, and the transmission across the sample scales as L~ x . 

In the quantized model (with % = 1) n takes integer values and the one-step evolution 
is specified by the unitary propagator: 

W = TE/o = [£ e- m2T ' 2 \n){n\+ £ \n){n\]e~ iV ^ . (2) 

nEM neZ\M 

The interacting dynamics defined by (2) and the free dynamics given by the propagator 
Uo = e~' v ( 6 } give rise to a well-posed scattering problem : the M0ller wave operators and 
the Scattering operator are mathematically well defined and so is the Scattering Matrix 
S a ,f3(ty a t fixed quasi-energy A. The scattering channels labelled by a are defined by 
free waves with wavenumbers 9 a satisfying kcos9 = A (moc?27r). For given A there are 
~ 2k l V channels, and associated with each channel 9 a there is a reversed channel 9 a * 
that corresponds to a wave propagating in the opposite direction. The Scattering Matrix 
is obtained in the standard way from distorted plane waves u+ a the latter are in turn 
solutions of the Lippman-Schwinger equation 

[1 - e a (U - e^+fV 1 " 1)]"+" = "o'" (3) 

where Uq'" are free waves. The Scattering matrix ^^(A) for a given quasi-energy A can be 
numerically computed with a great accuracy (unitarity error ~ 10~ 6 for sample lengths up 
to ~ 600). In our numerical computations we actually take V(9) = 2garctan(^ cos 9) that 
becomes kcos9 in the limit q — > oo,£ — > 0, 2q£ = k . The rank of the scattering matrix 
is ~ 2k I~k and becomes large in the classical limit (k — > oo,r — > 0, kr = const., L/k = 
const.). 
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Figure 1. Averaged squared moduli of S-matrix elements versus rescaled 
sample length from an ensemble of 100 matrices of rank N — 30, with 
JOT = 10, k = 45, L = 880, 30 < L < 240. Ope™ circles: 
diagonal reflection matrix elements. Open squares: nondiagonal reflection 
elements. Open triangles: transmission elements. Full symbols at scaled 
length ~ 4 are obtained by averaging over a quasi-energy interval. 

Characteristic length scales are given by the mfp ~ k and by the localization length I. 
The latter is semiclassically on the order of the classical diffusion coefficient D, which is 
in turn proportional to k 2 . The mfp is the one classically relevant scale and the classical 
transmission coefficient depends on the scaled sample length C = 2kL/irD ~ L/k, that 
is, on length measured in units of the mfp. The ballistic regime is defined by C ~ 1 and 
the classical diffusive regime by C » 1, meaning that a large number of collisions occurs 
inside the sample. In addition to this, the quantum diffusive regime requires L < I, and 
these two conditions can be met the more easily, the larger k. 

An ensemble of S-matrices is generated by shifting the sample along the n— axis. 
Under a well-known assimilation of the quantum standard map to a solid-state model of 
the tight-binding type, moving the sample to different positions is equivalent to taking 
different realizations of the disorder, while keeping the sample in a fixed position. In 
other models for chaotic scattering there was no "disorder averaging" and the statistics 
were generated instead by slightly changing the value of the energy. We have made a 
systematic use of ensemble averaging, but in a few test cases averages over small quasi- 
energy intervals were also computed and a close agreement with ensemble averages was 
found (compare for example the full symbols with the open ones in Fig.l). 
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The S-matrix has one symmetry expressed by the reciprocity relation S a p = Sp* a * 
Thanks to this property the matrix S a p = S a p* is unitary and symmetric; this very 
matrix was the object of our statistical analysis. In the rest of this Letter it will be 
denoted just S and no further reference to the original matrix will be made. In Ref.[7] 
one important prediction of the theory of irregular scattering was tested, concerning the 
correlation of S-matrix elements at different quasi-energies spaced by e; it was found 
that in the ballistic regime the squared modulus of such correlations depends on e in an 
approximately Lorentzian way, as expected on semi classical grounds [2]. 

In this paper we focus on the statistics of S-matrix at fixed quasi-energy. Fig. 1 
shows the average of squared moduli of three different classes of matrix elements : (a) 
diagonal elements (associated with reflection with unchanged velocity : open circles), (b) 
non diagonal elements (open squares) and (c) transmission elements (open triangles). For 
each class the recorded average was obtained by a double averaging process, first over all 
the matrix elements in the given class, and then over the matrix ensemble. The results are 
shown for values of the scaled length that, as shown in ref[7], correspond to the ballistic 
and the quantum ohmic regimes. 

COE is certainly not a proper reference ensemble for the diffusive regime, because 
transmission and reflection elements have systematically different magnitudes and thus 
violate the COE invariance properties. On the other hand, even in the ballistic regime fast 
direct reactions are present, connected with classical orbits that do not undergo significant 
randomization inside the interaction region. These direct reactions give a non-fluctuating 
contribution to the S-matrix and have to be filtered out somehow. 

To this end we used the quite general method described in ref.[5]: first we computed 
an ensemble-averaged (nonunitary) scattering matrix S = (5) that could be turned into a 
diagonal matrix S' by a suitable transformation, (S' = USU T , with U unitary and U T its 
transposed matrix). On applying the same transformation to the original matrix S a new 
symmetric unitary matrix 5" was obtained, which is diagonal "on the average". Finally, we 
investigated the statistics of the "fluctuation" matrix S? 1 = S' — S', searching for universal 
characters that could be modelled by Random Matrix Theory. In particular, we have 
compared our numerical results with the prediction of a general random matrix model 
which was found to satisfactorily reproduce the fluctuations observed in a dynamical 
model[5]. In the approach of Ref.[5] the average S-matrix is used as an input information 
in the building of a theoretical ensemble. The latter is uniquely defined by the coefficients 
T a = 1 — l^'ac-l 2 ; in particular, given the T a 's, theoretical values for the 2nd moments 
< 1 5^ | 2 > can be found by numerically computing certain multifold integrals [9]. In 
a number of ballistic cases (C ~ 1) we found the empirical distributions of S a p to be 
Gaussian with good accuracy. Empirical values of < \S^\ 2 > from an ensemble of 
100 matrices were found to agree with theoretical ones within an average 20% accuracy, 
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Figure 2. VaWaTice o/ diagonal elements of the fluctuating part of the 
S — matrix, for the same data as Fig.l and L=28. The labeling of matrix 
elements given on the horizontal axis is taken so that the coefficients T a 
monotonically increase from left to right. In order to facilitate comparison, 
a continuous line has been drawn through the theoretical data. 



both for diagonal and for off-diagonal elements. A comparison for diagonal elements 
is shown in Fig. 2. Although there is some qualitative agreement (in the non-diagonal 
case, the theoretical and the empirical 2nd moments remain within 20% of each other 
while changing over 15 orders of magnitude) the observed deviation is significantly larger 
than statistically expected; also, Fig. 2 gives a clear indication for the theoretical values 
systematically underestimating the empirical ones. In the diffusive regime the agreement 
becomes worse and worse, the discrepancies being particularly large at small values of 
the coefficients T a , where theoretical data decrease and empirical ones increase instead. 
However in the diffusive case the random model is not even theoretically adequate (it 
assumes all the bound states to be equally coupled, in a statistical sense). 

If the scaled length is increased still further, the localized regime is entered, where 
transmission elements become exponentially small. The statistics of the transmission 
coefficient T[12] in the localized regime presents some elements of interest. According 
to our data, not reported here, the distribution of T in this regime is approximately 
lognormal, consistently with theoretical predictions for quasi - Id Anderson insulators. 

In the far localized regime the S-matrix of rank 2N is approximately decomposed in 
two submatrices of rank N describing reflection on the left and on the right of the sample. 
In the limit of an infinitely long sample, these matrices become unitary and symmetric 
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Figure 3. Normalized distributions of the squared moduli of the S-matrix 
diagonal reflection elements (a) and of nondiagonal reflection elements (b)for 
a localized case, from an ensemble of 200 matrices of rank 8 and for 
kT — 5, k — 12.4, I ~ 50, L — 400, The s uperimposed 
theoretical distributions, as given in ref.flO], are 3/2\/l — X for case (a) 

and 3/2(1 - x) 2 2 F 1 (3/2, 1, 3, 1 - x) for case (b), where F is the 

hypergeometric function. 



themselves. In this regime we have found the statistics of the reflection submatrices 
to closely agree with COE predictions. In Fig. 3 the distributions of squared moduli of 
diagonal (a) and nondiagonal (b) reflection matrix elements are shown for a localized case 
i/L ~ 1/8, together with the theoretical distribution derived in ref.flO] (the latter is not 
an exponential one because of the small rank of the matrix, and in the nondiagonal case 
it involves an hypergeometric function). The average squared modulus of non diagonal 
elements is ~ 1/(N + 1) and that of diagonal elements is ~ 2/(N + 1). The factor 
2 accounts for backscattering[ll]. Here again the distributions were generated by the 
double averaging process described above. An important remark is that this agreement is 
here found in a highly nonclassical regime, dominated by the purely quantal localization 
effect. Indeed, one important prediction about S-matrix fluctuations in chaotic scattering, 
that was based on semi classical considerations [2], is here violated: the squared- modulus 
of the correlation of S— matrix elements at different quasi-energies is not of Lorentzian 
type. In Fig. 4 we show a typical correlation, 

Ca0[e) = (Mow (4) 
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Figure 4. Squared modulus of the normalized autocorrelation of a reflection 
matrix element, versus quasi-energy and for the same data as for Fig. 3. 

as a function of the quasi-energy e, for the same data as in Fig. 3. This dependence is 
clearly not a Lorentzian and is suggestive of a discontinuous derivative at e = in the 
limit of an infinitely long sample. We have no theoretical argument for such a behaviour. 
Discontinuous derivatives of the correlation curve at e = have been reported and ex- 
plained on semiclassical grounds in certain problems of chaotic scattering in the presence 
of classical nonexponential decay [13]. However semiclassical arguments not accounting 
for localization cannot work in the present case. 

In conclusion, we have numerically analyzed the fluctuations of the scattering matrix 
in a model which displays both classical chaotic diffusion and quantum localization. This 
model has been shown elsewhere to exhibit some important features of transport in dis- 
ordered 1-d solids; therefore the indications that emerge from our data have a bearing 
not only on the general problem of quantum scattering in the presence of classical chaos 
but also on the problem of transmission fluctuations. According to the results illustrated 
in this Letter, a rough qualitative agreement with the GOE model is observed in the bal- 
listic regime but systematic deviations exist. COE statistics emerge instead in the limit 
of extreme localization, where the problem becomes essentially one of reflection of waves 
by a disordered medium. In that case, however, the correlation of S — matrix elements 
at different quasi-energies is not a Lorentzian but it appears exhibit an angular point at 
the origin. It may be interesting to speculate whether this behaviour may be related to 
long time diffusive decay, but the presence of localization calls for great caution in that 
matter. 
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